Quantized Vortex States of Strongly Interacting Bosons in a Rotating Optical Lattice 
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Bose gases in rotating optical lattices combine two important topics in quantum physics: super- 
fluid rotation and strong correlations. In this paper, we examine square two-dimensional systems at 
zero temperature comprised of strongly repulsive bosons with filling factors of less than one atom 
per lattice site. The entry of vortices into the system is characterized by jumps of 2ir in the phase 
winding of the condensate wavefunction. A lattice of size L x L can have at most L — 1 quantized 
vortices in the lowest Bloch band. In contrast to homogeneous systems, angular momentum is not a 
good quantum number since the continuous rotational symmetry is broken by the lattice. Instead, 
a quasi-angular momentum captures the discrete rotational symmetry of the system. Energy level 
crossings indicative of quantum phase transitions are observed when the quasi-angular momentum 
of the ground-state changes. 



I. INTRODUCTION 

Rotating Bosc-Einstein condensates in dilute alkali 
gases provide a rich playground for the study of quantized 
vortices in superfluid systems. One major advantage is 
the ability to directly image vortex cores and study their 
static and dynamic properties. Vortices, first made by 
quantum state engineering of condensate wave functions 
[l|, |2j , are typically now produced by mechanical stirring 
of ultracold atomic clouds 0, 0] . Remarkable images 
of large vortex lattices containing more than one hun- 
dred vortices in an Abrikosov type triangular configura- 
tion have produced striking evidence for the superfluidity 
of Bose-Einstein condensed alkali gases 0, . 

Employing the analogy between the Hamiltonian for a 
two-dimensional electron gas in a strong magnetic field 
and that for a rotating atomic gas, it has been pointed 
out that the physics of the fractional quantum Hall effect 
(FQHE) should emerge when the number of vortices and 
the number of atoms become comparable 0, El E| • 
Achieving this regime experimentally is a significant goal 
of the field at this time. One of the reasons for this is the 
connection with many problems of interest in condensed 
matter systems where strongly correlated electron effects 
have been discussed and studied. The direct approach 
of spinning-up a Bose-Einstein condensate to reach the 
regime of strongly correlated effects is difficult because 
of the need to reach a parameter regime of low particle 
number per vortex and extremely low temperature 
for this reason the FQHE regime has yet to be achieved 
with cold atoms. 

Ultracold gases are typically dilute, and the interaction 
effects can be well incorporated by perturbation theory. 
However, it is possible to manipulate and enhance the in- 
teraction effects in a number of ways so that the pcrtur- 
bative treatment fails. One possibility is to increase the 
two-body scattering length via a Feshbach resonance. An 
alternative method is to modify the effective interactions 
through application of an optical lattice. An optical lat- 
tice is formed from an off-resonant light intensity pattern 
created by the interference of several laser beams. The 



atoms feel a potential proportional to the intensity of 
the light field. As the laser fields are made more intense, 
the interactions are enhanced because the atoms become 
more strongly confined in the lattice wells. This tech- 
nique was exploited in the theoretical and experimental 
study of the Mott-insulator to superfluid quantum phase 
transition in a non- rotating system [T2I . Theoretical 
studies have been done linking the Hofstadtcr butterfly 
[3 and the FQHE [H G3 with bosons in an optical 
lattice in the presence of an effective magnetic field. 

It is thus a natural question to pose as to whether one 
can combine the intriguing physics of the rotating gas 
with the enhanced interactions in an optical lattice. The 
motivation would be to move the regime of quantum Hall 
physics towards a parameter space that is experimentally 
achievable. In a first experiment with a rotating opti- 
cal lattice, Tung et al. [17| have recently demonstrated 
vortex pinning in a weakly interacting BEC. This was 
realized by passing a laser beam through a mask that 
contains holes arranged in a particular configuration and 
then focusing the laser beams to form the lattice inter- 
ference pattern. The two dimensional optical lattice is 
rotated by spinning the mask. 

In this paper, we formulate a theoretical description 
of bosons in rotating optical lattices and illuminate the 
connections with vortex physics, building upon results 
presented in an earlier article [18j]- Certain aspects of 
this problem such as single vortex formation |l9l| and vor- 
tex pinning [2(j have been theoretically explored for high 
filling factors (large number of atoms per site). Burkov 
et. al. [2]J have shown the formation of delocalized vortex 
clusters in lattices containing superfluids in the presence 
of an effective magnetic field created by modulating the 
optical lattice. We consider strongly repulsive bosons in 
a small 2D rotating square optical lattice with filling fac- 
tors less than unity. In this regime, exact solutions are 
tractable and indicate key properties of larger systems. 

The system can be studied using a modified Bose- 
Hubbard Hamiltonian with a complex, site-dependent 
hopping term which is sufficient to describe physics in 
the lowest Bloch band. The lattice breaks the continuous 
rotational symmetry associated with the angular momcn- 
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turn operator. Instead, it possesses a discrete rotational 
symmetry, and the generator of this discrete rotation, the 
quasi-angular momentum, plays an important role. The 
square lattice is four-fold rotationally symmetric and the 
quasi-angular momentum operator generates rotations in 
steps of 7r/2. 

The entry of vortices into the system is marked by 2tt 
jumps in the phase winding of the condensate wavefunc- 
tion. We show that a maximum vorticity of 27r(L — 1) 
is possible for the lowest band in a lattice of size L x L. 
Changes in the quasi-angular momentum of the ground 
state are associated with energy level-crossings as a func- 
tion of changing angular velocity. For filling commensu- 
rate with the symmetry of the system, the quasi-angular 
momentum is zero at all angular velocities. In this case, 
avoided energy level crossings of the ground-state are ob- 
served as a function of lattice angular velocity. For in- 
commensurate filling, there are energy level crossings for 
many particles at zero temperature. Since these corre- 
spond to a symmetry change in the ground state as a 
function of a Hamiltonian parameter at zero tempera- 
ture, and the property holds for systems of arbitrary size 
and particle num ber . these critical points are quantum 
phase transitions |27j |. 

This paper is structured as follows. Section [H] pro- 
vides the theoretical framework by sketching the deriva- 
tion of a modified Bosc-Hubbard Hamiltonian for describ- 
ing lowest-band physics and discusses its regime of va- 
lidity, along with the principal observables. Section IIIII 
describes general characteristics of the system such as 
the Mott-insulator/superfluid phase diagram and the ef- 
fect of strong repulsive interactions. The behavior of 
the density distribution as a function of rotation is il- 
lustrated. Section II VI presents a characterization of the 
quasi-angular momentum. Sections Ivl and IVll discuss re- 
sults for a single particle and many particles in the sys- 
tem, respectively. Section IVIII summarizes the main re- 
sults of this paper. 



II. MODIFIED BOSE-HUBBARD 
HAMILTONIAN FOR ROTATING LATTICES 
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FIG. 1: Schematic for a 4 x 4 lattice rotating in the counter- 
clockwise direction. Sites labeled i and j are nearest neigh- 
bors. The rotation-driven hopping between neighbouring sites 
is governed by the parameter Kij = /3(xiVj — XjUi)/d where 
j3 is given by the overlap integral Eq. © and the second fac- 
tor is the perpendicular distance of the line joining the two 
neighboring sites from the center of rotation. 



The Hamiltonian in the rotating frame reads 
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where M is the single particle mass and g is the coupling 
constant for repulsive two-body scattering in a dilute gas. 
The bosonic field operator $ obeys the commutation re- 
lationship 



<J>(x),$(x') 



<5(x 



). Particles can be 

described using an orthonormal Wannier basis (x) . 
Here, i indexes the N sites on the lattice and p denotes 
the band index ^3] • If the energy due to interaction and 
rotation is small compared to the energy separation be- 
tween the lowest and first excited band, the particles are 
confined to the lowest Wannier orbital. We shall consider 
this regime only and henceforth will drop the band index 
p. The field operator <E> can be expanded in terms of this 
Wannier basis, Wj(x), and the corresponding site-specific 
annihilation operators, CLi , as 



A. Hamiltonian 

The energy for a fixed number of bosons in an optical 
lattice can be broken down into three components corre- 
sponding to the kinetic energy, the potential energy due 
to the lattice, and the interaction energy between bosons. 
Using a standard procedure , the Hamiltonian in the 
reference frame rotating with angular velocity f2 about 
the z-axis is H = Hq — j dx^flL z ^, where Ho is the 
Hamiltonian in the laboratory frame and L z is the an- 
gular momentum. This coordinate transformation facili- 
tates the calculation of the ground state in the laboratory 
frame since it renders the Hamiltonian time-independent. 



N 

$(x)=J>Wi(x). (2) 

i=i 

Alternatively, a rotation dependent phase can be ascribed 
to each Wannier basis element and the field operator <3E> 
expanded accordingly. Comparisons with this approach 
are made in the next subsection. 

In the tight binding regime, tunneling between sites 
which arc not nearest neighbors can be neglected. The 
interaction between particles on nearest neighbor sites 
can also be neglected. Using this approximation and sub- 
stituting Eq. (J2J) into Eq. yields the modified Bose- 
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Hubbard Hamiltonian 

H = - 



(a\a,j + aiafj + e^hi + ^ ^ - 1) 
-iTtil^Kij {alhj - ^afj , (3) 

(i,3) 

where i and j are site indices, indicates that the 

sum is over nearest neighbors, and hi is the number op- 
erator for site i. The first three terms are common to 
the well-studied Bose-Hubbard Hamiltonian for particles 
in a stationary lattice [23j. The parameters t and e are 
integrals describing hopping and onsite zero-point energy 
respectively: 



t EE / dxWT(x) 

/ dxW*{x) 



-— v 2 + y^(x) 

2M 1 ' 

2M y ' 



W 3 -(x), (4) 
W,(x). (5) 



Wannier functions along the a; and y directions can be de- 
coupled for a square lattice, and accordingly, the integrals 
become one-dimensional. The third term in the Hamilto- 
nian, Eq. J2J), describes the interaction between particles 
on the same site. For an s-wave scattering length a s 

MM, 



u 



4ira s h 
M 



J dxlVt^x)! 4 



(6) 



The last term in Eq. © is the modification due to the 
rotation and favors hopping along one azimuthal direc- 
tion. Kij is a product of the azimuthal overlap integral, 
/3, which is dependent on the geometry and form of the 
lattice, and the perpendicular distance of the line joining 
sites i and j from the center of rotation (Fig. [IJ, 



K . P( _ X 



(7) 



Here, (xi,yi) are the coordinates of the i th site with the 
origin located at the center of rotation. The lattice spac- 
ing is d and 







dxW*(x-d)d x Wj(x) 



(8) 



We numerically calculate t and j3 for a sinusoidal lattice 
potential, V = Vb(sin 2 (7ra;/d) + sin 2 (it y / d)) using Math- 
ieu functions. The results are plotted in Fig.|2|as a func- 
tion of the lattice depth up to a very tight confinement of 
Vq/Er = 20, with E R = h 2 TT 2 /2md 2 denoting the recoil 
energy. 

The eigenstates of the Hamiltonian can be written in 
the form 



l*> = 



X! c {n,}\ni 1 n 2l 

{»«} 



,n N ) . 



(9) 
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FIG. 2: Overlap integrals t (Eq. gj) and (Eq. ®) 
for a standing wave optical lattice described by V = 
Vo(sin 2 (jrx/d) + sin 2 (ivy /d)). The lattice depth is given in 
units of the recoil energy Er = Tl 2 tt 2 /2md 2 . (a) The hopping 
parameter t decreases exponentially as a function of Vq. (b) 
pdER/t as a function of lattice depth. In the tight-binding 
regime, j3d scales roughly linearly with t/En. In this paper 
we use j3 — 4.93t/ End and t — 0.02Er , corresponding to a 
standing wave optical lattice of depth Vq/Er = 10. 



Here indicates the total number of sites and {n^} indi- 
cates the set of all possible products of number states con- 
strained by the total number of particles, i.e. rii = n. 
In this paper, we use the truncated set of Fock states 
1 1 0), to describe the number of particles at each 
site |l6t l25j . This corresponds to assuming a regime 
of strong repulsive interactions, i.e., that of hard-core 
bosons. Note that in a ID lattice this approach is equiv- 
alent to mapping bosonic operators onto fermionic ones 
via the Jordan- Wigner transformation j2Q, H3| • A test of 
the regime of validity of this approximation is provided 
in Section fill Al The Hamiltonian is constructed using 
this basis and diagonalized to find the ground state en- 
ergy eigenvalue and eigenstate. A set of tools, which are 
described in Section III CI are then used to analyze the 
ground state. Note that solving for the Hamiltonian in 
Eq. © is equivalent to putting the 2D lattice inside a 
box with infinite potential walls. This effect leads to a 
number density distribution that is peaked in the center 
for zero rotation. 



B. Regime of Validity of the Hamiltonian 

Many aspects of studying the characteristics of bosons 
in a rotating optical lattice map onto the extensively 
studied problem of Bloch electrons in the presence of a 
magnetic field [2^, l29l I30L l3lj . This subsection makes a 
connection with the electron problem, while exploring, in 
parallel, the limitations of using the lowest band Bose- 
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FIG. 3: Projection of the ground state obtained using imagi- 
nary time propagation (ITP) onto the Hilbert space spanned 
by the eigenvectors obtained using H (solid line) (Eq. ©) and 
H2 (dashed line) (Eq. JlH l for a 2 x 2 lattice with a lattice 
depth of Vo = WEr. f2 is in units of the recoil energy. The 
overlap is good even up to Q = Er/H ~ 50t where t is the 
hopping energy. 



Hubbard Hamiltonian. 

The Hamiltonian described in Eq. only takes into 
account Wannier orbitals from the lowest band — an ex- 
cellent approximation for small CI and a deep lattice. An 
alternative formulation can be obtained by using Eq. Q 
in conjunction with a different Wannier basis given by 



W/(x)=exp 



iM 



f A(x') • dx' 



Wi(x), (10) 



resulting, for the case of a single particle in which interac- 
tion effects arc not present, in the following Hamiltonian, 



Ho 



-E 

(id) 
x exp 



1 



t + -mflH 



2 " * ojo,- 



A(x') • dx' 



h.c. 



FIG. 4: Spatial ground state number density (left) and phase 
information (right) for one particle in a 2 x 2 sinusoidal lattice 
with Q — 0.5Er/Ti and Vo = 10-E.r obtained using (a) Imagi- 
nary time propagation, (b) Hamiltonian Hi (Eq. (J3J) and (c) 
Hamiltonian H 2 (Eq. Illl l. 



E 



e - -mn 2 rf - -mn 2 e' 



(ii) 



Here, i' and e' refer to integrals similar to those defined 
for t and e in Eqs. |J3J and (J5J. A(x) = f2 x x is the ana- 
log of the magnetic vector potential. For more than one 
particle, this Hamiltonian must be extended to include 
the effects of interactions. Note that H2 and W- (x) have 
forms similar to those used traditionally in the treatment 
of Bloch electrons in the presence of magnetic fields. A 
similar formulation has also been used to study bosons 
in an optical lattice in the presence of an effective mag- 
netic field [Til ITU and for bosons in a rotating optical 
lattice 0. 

In this subsection, the two approaches are compared 
with results for a single particle obtained from imagi- 
nary time propagation (ITP) for one particle in a 2 x 2 
lattice. This analysis gives rise to three findings: (1) The 
ground state of the Hamiltonian described in Eq. (J3J) no 
longer depends on the increase in f2 once a maximum 
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phase difference of 7r/2 between neighboring sites has 
been reached, i.e., all the vortex entry transitions possible 
within the lowest band have occurred. For one particle in 
a 2 x 2 lattice, the corresponding maximum phase wind- 
ing is 2tt. This limitation does not apply to the Hamil- 
tonian Eq. (2) The Hilbert spaces spanned by the 

cigcnstatcs of both Hamiltonians Eq. © and Eq. 
capture most of the exact ground state wavcfunction for 
hil < Er. Note that TiVt ~ Er is large from an ex- 
perimental point of view. For the case of one particle 
in a 2 x 2 lattice, the projection of the exact wavefunc- 
tion on either Hilbert space is > 90% for Ml ~ Q.5Er 
(Fig. Both approaches yield accurate density pro- 
files for large Q (~ Er/K) but differ from the ITP-result, 
and from each other, with regard to the velocity pat- 
tern. Note that H and H2 involve different approxima- 
tions to the phase gradient. The Hamiltonian H allows 
for phase changes only in the region of overlap of next- 
neighbor Wannier functions Wi(x), i.e., yields a uniform 
phase around the site center (Fig. Efb)). The Hamilto- 
nian H2 requires phase gradients to be proportional to 
f2 and allows for non-zero phase gradients within each 
well (Fig. U c ))- (3) The lattice rotation frequencies at 
which the first vortices appear, as will be discussed later, 
are slightly different in the two cases over the range of 
interest due to the different influence of higher bands in 
the three formulations. The remainder of our work is pri- 
marily concerned with the states of the system at low O 
and since these are well captured by the simpler lowest 
band Hamiltonian described in Eq. J2J, we use it for the 
rest of the paper. 

C. Toolkit 

We evaluate six quantities to characterize the behav- 
ior of the system: (1) energy, (2) site number density, (3) 
intcrsite current in the rotating frame, (4) average angu- 
lar momentum, (5) quasi-angular momentum eigenvalues 
and (6) phase winding of the condensate wavcfunction. 

The ground state energy is obtained as the lowest 
eigenvalue of the Hamiltonian. The site number den- 
sity is the expectation value of the number operator: 
i%i = (a\a,i). The expectation value for the current 
flowing from site i to a neighboring site j in the rotating 
frame is obtained using the continuity equation, 

= - ~ toKijfeaj + ajaj), (12) 

where the current is in units of t/h. Hij in Eq. 1)1 2[l is the 
part of the Hamiltonian operator relevant to sites i and 
j. Since the number density on any site i is constant for 
any steady state solution, the algebraic sum of currents 
associated with any site i is zero. 

The derivative of the energy with respect to the angu- 
lar velocity Q — keeping all other Hamiltonian parame- 



ters constant — gives direct access to the average angular 
momentum, 

BE 

(l.) = -35. (13) 

Quasi-angular momentum is calculated based on the 
four-fold rotational symmetry in the Hamiltonian created 
by the square lattice, i.e., rotating the system by 90 de- 
grees has no effect on the Hamiltonian. This is discussed 
in detail in Section llVI Here, we give a brief description 
of the way we calculate the quasi-angular momentum. 
A discrete rotational symmetry operator, R(7r/2), can 
be constructed such that R(7r/2) acting on a wavcfunc- 
tion rotates that wavefunction by ninety degrees. R(7r/2) 
commutes with the Hamiltonian and therefore shares si- 
multaneous eigenfunctions. Applying this operator four 
times corresponds to a rotation of 27r, bringing the wave- 
function back to its original state: 

R* = 

R% = $^r = e'" ll/2 , me {0,1, 2, 3}. (14) 

Hence, the eigenvalues of R(ir/2) arc defined by the 
quantized dimensionless quasi-angular momentum m. 
From an operational standpoint, we apply a rotation 
R(7r/2) to \& and read out the eigenvalue, r. Note that 
the discrete rotational symmetry operator can be gener- 
alized to a n— fold rotationally symmetric system. 

The condensate wavefunction is the eigenfunction cor- 
responding to a macroscopically large eigenvalue of the 
one-body density matrix G^, with all other eigenvalues 
being non-macroscopic [32,|33|]. In usual tensor notation, 
G-j' = (&jO,i). The phase of the condensate wavefunc- 
tion describes the superfluid properties of the system. In 
the small systems we consider in this paper, this is not 
a rigorous definition. However, a meaningful condensate 
wavefunction can still be obtained in this way since in the 
superfluid regime one of the eigenvalues of the one-body 
density is always significantly larger than all the others, 
even for a very small number of particles. The phase 
winding around the perimeter of the condensate wave- 
function, O c f, when divided by 2-7T gives the vorticity of 
the system. The subscript cf indicates that the phase 
winding refers to that of the condensate wavcfunction. 

III. GENERAL CHARACTERISTICS 

This section describes four general features of bosons 
in rotating lattices: (1) interaction effects, (2) the Mott 
insulator/superfluid quantum phase diagram for the sys- 
tem, (3) number density depletion at the center for odd 
lattices, and (4) number density distribution for even lat- 
tices. The first subsection lays out the justification for 
our choice of studying lattices using a truncated Fock 
space as mentioned above in connection with Eq. 10 ■ 
The second subsection connects with existing under- 
standing of the phase diagram for the Bosc-Hubbard 
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FIG. 5: Current J12 between adjacent sites of a 2 x 2 unit cell 
lattice as a function of total number of particles n for different 
values of the ratio between repulsive interaction and hopping 
U/t at large rotation (HQ ~ 0.6Er S> t). For weak interac- 
tion, U = 0.5t, the particles can freely cross each other and 
the current is proportional to the number of particles. As the 
interaction increases, the current per particle drops as a func- 
tion of filling. For U = 10004, the current for three particles 
is the same as that for one particle, indicating a particle-hole 
symmetry. The current was evaluated using a four-state Fock 
basis on each site, which allowed three particles per site. 



model. The third subsection discusses our choice of lat- 
tices with an even number of sites and the last describes 
number density rearrangement with vortex entry. 



A. Effect of interaction 

Interaction between bosons inhibits current flow by 
making it difficult for particles to cross each other. To 
demonstrate this we consider currents in a 2 x 2 unit 
lattice using a four-state Fock basis on each site. Fig- 
ure is a plot of the current in the rotating frame be- 
tween two neighboring sites as a function of filling for 
different interaction strengths at fixed angular velocity 
(hfl ~ 0.6-Er 3> t). For weak interactions (U = 0.5i), 
the current is proportional to the number of particles 
as they can flow independently of each other. However, 
the current per particle drops with increasing interaction 
and filling. At large interaction (U = lOOOi), the cur- 
rent for three particles (one hole) is the same as for one 
particle in the system. This particle-hole symmetry is 
characteristic of the regime where bosons are impenetra- 
ble, i.e., of the regime where the two-state approximation 
applies In fact, currents calculated for U > lOOt 

using the two-state approximation coincide with those 
obtained with a larger Fock space. The main results of 
this paper are obtained assuming the atoms to be im- 
penetrable and hence are expected to be quantitatively 
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FIG. 6: (a) Mott-Insulator zone boundaries for bosons in a 
2x2 lattice as a function of chemical potential fi and of the 
real and imaginary parts of the complex hopping parameter 
described in Eq. 1161 . The volume inside the shaded surfaces 
corresponds to Mott-insulating states with average fillings of 
one and two particles per site in the lower and upper volume 
respectively, (b) The x — z cross-section of the phase diagram 
is the same as that obtained using the Bose-Hubbard Hamil- 
tonian for a non-rotating system, (c) The y — z cross section 
of the phase diagram is similar to the x — z cross section. Note 
that the shape of the Mott zones will change on going beyond 
the two-state approximation 13^1 . 



accurate in the regime U > lOOt for fillings < 1. 



B. Phase diagram 

The phase diagram for the non-rotating Bose-Hubbard 
model is obtained in the grand-canonical ensemble by 
adding a term —fx hi to the Hamiltonian. The phase 
diagram separates into two regions — the Mott insulator 
(MI) lobes with commensurate filling (integer number 
of atoms per site) and the supcrfluid (SF) regions with 
incommensurate filling. The phase diagram for the non- 
rotating lattice was first studied in Ref . j2^| ■ 

The introduction of rotation makes the hopping energy 
parameter t complex and site-dependent as can be seen 
by rewriting Eq. @ in the grand canonical description 
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y/d x /d 



FIG. 7: Number density distribution for one particle in a 3 x 3 
lattice of depth Vb = WEn at a lattice rotation of h£l — Er. 
There is almost complete number depletion in the central site 
because it coincides with the center of rotation and hence the 
vortex core. 



after combining the first and last terms. 



(a) (c) 




y/d -2 -2 x / d y/d -2 -2 M 



FIG. 8: Number density distribution for one particle in a 4 x 4 
lattice at Vb = WEr. (a) TiQ. = O.IEr - rotation has yet to 
enter the system. The center-peaked distribution is due to 
the infinite potential walls at the perimeter of the lattice, (b) 
frtt = 0.2E R (c) TiQ. = QAE R (d) hQ, = 0.8E R . At large 
rotation, particles get pushed to the outermost sites creating 
a washing machine effect. 



H = (tijolaj + l]n,n\) + (e - fx) hi 

+ -J2Mni-l), (15) 

i 

Uj = t + ihnK i3 . (16) 



The surface in Fig. represents the boundary between 
the Mott insulator (MI) and superfluid (SF) regions as 
a function of the chemical potential and of the real and 
imaginary parts of the hopping parameter ty. The zero 
point energy e in each well has been set to zero since 
it gives rise to an irrelevant overall shift. Note that 
the shape of the tip of the Mott zones in Fig. differs 
slightly if more than two Fock states are used to repre- 
sent each state Q. At fl = (Im(%) = 0, Fig. H^b)), 
the phase boundary matches that for the non-rotating 
Bose-Hubbard model in the two-state approximation. A 
similar diagram is obtained when setting Re(Ay) = 
(Fig-I3 c ))- In general, this similarity need not be present 
for larger systems where Uj is not the same for all pairs 
of nearest neighbors as in the example of a 2 x 2 lattice. 
Note that even though Im(t^ ) can be varied freely via f2 
for a given t and non-zero Kij, the reverse is not true. 
This is because for a particular realization of a lattice, the 
two overlap integrals t (Eq. (QJ) and (3 (Eq. ©) are both 
fixed by the lattice depth. Hence, a fixed Im(ty-) = flKij 
implies a fixed Ke(tij) = t. In the case of a standing wave 
optical lattice in the tight binding regime, Kij is approx- 
imately proportional to t/En (see Fig- EJ - This makes 
the plane Re(ty ) = of the phase diagram inaccessible. 



C. Differences between even and odd lattices 

In the study of rotating lattices, the position of the 
center of rotation gives rise to an important distinction 
between two kinds of lattices: lattices with an even num- 
ber of sites (e.g., 2 x 2, 4 x 4) and lattices with an odd 
number of sites (e.g., 3 x 3, 5 x 5). If assumed to be 
at the center of the system, the axis of rotation passes 
through a peak in the lattice potential in the case of even 
lattices while it passes through the central site for odd 
lattices. It is useful to briefly touch upon two interest- 
ing aspects of the ground state solution for odd lattices: 
(1) there is nearly complete number density depletion in 
the center site when rotation enters the system, as would 
be expected at the vortex core in a continuous system 
(Fig. 0; and (2) there are no currents along the radial 
direction (K^ = 0) so that particles only hop along the 
azimuthal direction. Even lattices have no sites which 
are nearest neighbors in a strictly radial direction, but 
the main results for particles in a rotating even lattice 
can be mapped onto those for particles in an odd lattice. 
Only the even lattice is discussed for the rest of the paper 
since it captures the most important physical features. 

D. Number density distribution 

The number density distribution is obtained by eval- 
uating the expectation of the site-specific number oper- 
ator hi. Figure [S] describes the number distribution for 
one particle in a 4 x 4 lattice with phase windings of 0, 
27r, 47r, and 6tt. For small angular velocities, there is no 
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FIG. 9: (a) Energy-quasi angular momentum dispersion rela- 
tionship for a 12-site ring. Positive and negative m\2 values 
indicate rotation in opposite directions. The velocity of the 
particles in any state is given by the slope of the dispersion 
curve at the point. Accordingly, the velocity is maximum at 
m = 3 (b) Energy — quasi-angular momentum dispersion re- 
lationship for a 12-site ring with a small four-fold periodic 
potential. Note the lines joining states are obtained by ex- 
trapolating for an infinite system. 



effect of rotation on the ground state of the system and 
the number density is center-peaked (Fig. Eta))- In the 
continuous limit of the lattice spacing becoming infinites- 
imal, this mirrors the number distribution for particles in 
a 2D box with infinite potential walls. The distribution 
changes each time a vortex enters the system. For large 
rotation (Fig.|SJi), we observe a washing machine effect, 
as the number density gets concentrated at the perime- 
ter. The inner sites have a small non-zero number density 
in this limit. 



IV. QUASI-ANGULAR MOMENTUM 

The presence of the lattice breaks the continuous ro- 
tational symmetry of the system. The eigenvalues of the 
angular momentum operator are therefore no longer good 
quantum numbers because the rotational symmetry as- 
sociated with L has been replaced with a discrete rota- 
tional symmetry. In this section, ideas of discrete transla- 
tional symmetry and Bloch's theorem are mapped onto a 
discrete rotational symmetry problem to generate quasi- 
angular momentum states. Exact results are presented 
for the modified Bose-Hubbard model in the context of a 
single particle one-dimensional ring, and connections are 
made with the square lattice. 

Consider a one-dimensional lattice with periodic 
boundary conditions, i.e., a ring lattice of N sites. A ro- 
tation o{2it/N leaves the system invariant and hence the 
rotation operator R (2w/N) commutes with the Hamilto- 
nian. This is also true for a square ring because the site 
dependent parameter in the Hamiltonian (Eq. @) 



depends on the perpendicular distance of the line con- 
necting two nearest neighbor sites from the center of ro- 
tation. The energy eigenstates can be labeled using the 
eigenvalues of R (2-k/N): 

R{2ir/N)\m) = e l27Tm/N \m) : to g {0, . . . , N- 1}, (17) 

where the exact eigenvectors shown can easily be derived 
by expanding |m) in the Fock basis and demanding peri- 
odic boundary conditions. 

At this point, it is useful to make a connection with 
conventional Bloch theory. R(tt/2) is analogous to the 
discrete translation operator TfcQ for a stationary one- 
dimensional lattice of period d |35j, 

T(d)V( Xl ,...,x n ) =e i * d V[x 1 ,...,x n ), (18) 

where ^(xi, . . . , x n ) is an eigenfunction of the translation 
operator. The eigenvalues of T(d) arc described by the 
quasi-momcntum q. In a way exactly analagous to that 
of quasi-momentum Bloch states for a discrete transla- 
tion operator, we can identify these m-values in Eq. H17|) 
as quasi-angular momenta. Note that this discussion so 
far is completely general and applies to both the single- 
particle and many-particle cases. 

To illustrate the role of the quasi-angular momentum 
and the connection with the quasi-momentum in systems 
with discrete translational symmetry, consider one par- 
ticle in a 12 site static ring. Each of the sites is indexed 
by an azimuthal coordinate, fa. The energy spectrum 
takes on the well-known dispersion relation observed for 
the lowest Bloch band of a particle in a ID lattice with 
periodic boundary conditions in the tight-binding regime 
(Fig. |SJa)). Since the system has 12-fold symmetry, the 
quasi-angular momentum, to, can take on 12 possible val- 
ues. The slope of the energy plot provides the velocity of 
the particle. As discussed before (Section [DJ, rotation is 
introduced by adding a term, —(hVL/i)d^, to the Hamil- 
tonian in order to obtain the ground state in rotating 
frame coordinates. As Q is ramped in a particular direc- 
tion, the ground state quasi-angular momentum changes 
from to = to to = 3 in steps of 1 (not shown here). 
The to = 3 state corresponds to the maximum slope of 
the dispersion and the largest particle velocity. A par- 
ticle in the m = 4 state has the same velocity as the 
m = 2 state but with a higher energy. Quasi-angular 
momenta to = 7, . . . , 11 correspond to to = —5, . . . , —1 
and describe circulation in the opposite direction. This 
is described by the C12 point symmetry group. 

Consider now a 12-site ring perturbed by a four- fold 
symmetric periodic potential. The Hamiltonian, H = 
H12 + V, is the sum of two terms, the 12-site Hamiltonian 
H12 which has a 12-fold rotation symmetry and a poten- 
tial V which has a four-fold rotation symmetry. Fig- 
ure E£b) is the energy dispersion relation as a function of 
quasi-angular momentum for the 12-site lattice ring with 
this small four-fold symmetric potential. Since the po- 
tential increases the rotational symmetry from d = 2n/12 
to d = 2tt/A, the Brillouin zone is narrowed down to 
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m = — 2, . . . , 2. Three energy bands are created in place 
of one. States on adjacent bands with the same m value 
— for example mn = —2 and m\2 = 2 — are mixed by 
the four-fold symmetric potential V, thereby leading to 
an energy gap at this m value. 

An analogous situation occurs when we try to qual- 
itatively understand the properties of a 4 x 4 lattice, 
which has 12 sites on the boundary. We adopt a per- 
turbativc approach by breaking the system into two non- 
interacting 12-site and 4-site rings and considering an 
interaction between them, i.e., H = H\2 + H4 + V. The 
interaction with the four-site ring breaks the 12-fold sym- 
metry of the outer ring, reducing it to a four-fold discrete 
rotational symmetry. 

The above example illustrates that a particle in a 
square lattice is characterized by a four-fold discrete 
rotational symmetry. The same symmetry considera- 
tions hold for many particles in the system. Hence, 
the many-body eigenstates are quasi-angular momentum 
states with m e {-2, 1, 0, 1, 2}. 

In the following we show how the rotation of the lat- 
tice leads to a change in quasi-angular momentum in the 
groundstate of the system in the single-particle case (Scc- 
tion^l and in the many-body case (Section lVl|) . In addi- 
tion, we show how these transitions affect other proper- 
ties of the system, such as its average angular momentum 
and its vorticity. 



V. ONE-PARTICLE ANALYSIS 

This section examines the response of one particle to 
lattice rotation. The advantage in first considering only 
one particle is that it allows one to distinguish general 
characteristics of the systems from effects due to interac- 
tion. 

Figure ^| describes the response of the system as a 
function of the angular velocity CI. Figures HOf al and 
rrUT b) show that for increasing Q the ground state energy 
in the rotating frame Eq decreases with discontinuous 
derivative as different states become energetically favor- 
able. Note that the fact that the eigenstates are not 
eigenstates of angular momentum is explicitly illustrated 
here since (L z ) takes on non-quantized values. 

The abrupt changes in average angular momentum are 
connected to changes in the quasi-angular momentum m 
of the groundstate as seen in Fig. llUf c) . Since the lattice 
has four-fold rotational symmetry, the values that the 
quasi-angular momenta can take on arc m 6 {0, 1, 2, 3}. 
For additional transitions, m repeats itself in behav- 
ior analogous to that of linear quasi-momentum as one 
crosses the first Brillioun zone. 

The changes in quasi-angular momentum are associ- 
ated with changes in the phase winding of the single par- 
ticle wavefunction. The phase winding jumps by 2tt 
each time the quasi-angular momentum of the ground- 
state changes (Fig. HUf d)). The maximum phase wind- 
ing of 147T for an 8 x 8 lattice corresponds to a max- 
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FIG. 10: One particle in an 8 x 8 lattice, (a) Ground state 
energy, Eo vs. fi. (b) Average angular momentum. (L z ) vs Q. 
Note that the expectation value of angular momentum (see 
Eq. 1131 ') is not quantized, (c) Quasi-angular momentum m 
vs. Q. In direct analogy to quasi-momentum values for linear 
lattices, m repeats itself, (d) The phase winding, Q/2n, vs. 
f2. For the 8x8 lattice, the maximum phase winding is 7 x 2n. 



imum phase difference of 7r/2 between any two sites on 
the lattice boundary. A difference of 7r/2 between two lat- 
tice sites corresponds to the condition for the maximum 
current attainable within the lowest band Bosc-Hubbard 
model. This result can be generalized to a lattice of size 
L x L. The number of sites on the circumference of the 
lattice is A{L — 1) and for a phase difference of tt/2 be- 
tween two adjacent perimeter sites, the maximum phase 
winding around the circumference is 2tt(L — 1) within the 
lowest band. Since the notions of an order parameter and 
of superfluidity do not apply to single particle systems, 
states with non-zero phase winding can not be referred to 
as quantized vortices. However, as will be shown in the 
following, these single particle results extend to many- 
particle systems in a straightforward manner. 
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FIG. 11: Quasi-angular momenta for one through four 
strongly repulsive bosons in a 6 x 6 lattice. For multiple par- 
ticles in the lattice and increasing rotation, the quasi-angular 
momenta cycles through values given by m = nl mod 4, I = 
0, 1, 2, 3, where n is the number of particles, (a) n = 1 : m — 
0,1,2,3,0,1 (b) n = 2 : m = 0,2,0,2,0,2 (c) n = 3 : m = 
0,3,2, 1,0,3 (d) n = 4 : m = 0. 



VI. MANY-PARTICLE ANALYSIS 

This section probes the effects of strongly repulsive in- 
teraction when the number of particles in the system 
is greater than one. The symmetry considerations dis- 
cussed in Section IIVI are first tested for many parti- 
cles. The effect of a symmetry-commensurate filling is 
explored by considering two different systems - (1) four 
particles in a 4 x 4 lattice and (2) five particles in a 4 x 4 
lattice. These two systems correspond to fillings com- 
mensurate and incommensurate with the four-fold sym- 
metry of the lattice, respectively. 

As discussed in Section IIVI the four-fold rotational 
symmetry allows labeling of the many-body states by 
their quasi-angular momenta. In a static lattice, the 
groundstate is always characterized by m = 0. This value 
may change when the lattice is rotated at angular veloc- 
ity Q, as has been demonstrated in the single-particle 
case. In contrast to the single-particle case, the many- 
body groundstate does not necessarily cycle through all 
possible m- values as fi is increased. Instead, the values 
of quasi-angular momenta it can take on depend on the 
number of particles. As illustrated in Fig.[nj we see that 
for n particles in the system, the quasi-angular momen- 
tum of the ground state cycles through values satisfying 
the relation 

m = nl mod 4 , (19) 
where I S {0,1,2,3}. The validity of this expression has 
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FIG. 12: Four particles in a 4 x 4 lattice, (a) Number of 
vortices, Q c f/2n, vs. fl. Three vortices enter the 4x4 lattice. 

(b) Quasi-angular momenta, m vs. Q. The symmetry of the 
ground state as indicated by the quasi-momentum m = 
does not change even with three vortices entering the system. 

(c) Zoom-in of lowest two energy levels around the entry of 
the second vortex shows an avoided energy level crossing. The 
mixing of states is possible because the ground states on either 
side have the same discrete rotational symmetry. 



been verified both numerically for various lattice sizes 
and particle numbers and analytically within a Jordan- 
Wigner transformation approach to hard-core bosons in 
a ring |36l |. Hence, only for odd n does the quasi-angular 
momentum of the groundstate cycle through all values of 
m. Cases in which the particle number is commensurate 
with the four-fold symmetry n = 4, 8, . . . are of particular 
interest in that these systems always stay in an m = 
state. Note that a simplistic explanation for the validity 
of Eq. l(T§)l is obtained if all particles occupy a conden- 
sate mode with quasi-angular momentum m = 0, 1, . . ., 
yielding a total quasi-angular momentum equal to when 
?i = 4. 



A. Symmetry-commensurate filling 

When the number of particles is commensurate with 
the four-fold rotational symmetry, the groundstate al- 
ways has zero quasi-angular momentum. This does not 
exclude the entry of quantized vortices into the system. 
To give an example, we analyze the case of four particles 
in a 4 x 4 lattice. The largest eigenvalue of the ground 
state one-body density matrix is found to be 60 — 74% 
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FIG. 13: Five particles in a 4 x 4 lattice, (a) Number of vor- 
tices 0C//2-7T vs. fi. As before, a maximum of three vortices 
enter the system (b) Quasi-angular momenta m vs. Q. Since 
the filling is incommensurate with the symmetry, the quasi- 
angular momenta takes on values m — 0, 1, 2, 3 (Eg. 1191 . (c) 
Zoom-in of lowest two energy levels around the entry of the 
second vortex shows an energy level crossing as the ground 
state symmetry changes as a function of a parameter in the 
Hamiltonian, fl. 



of the total particle number. Since all other eigenvalues 
are significantly smaller, this number is large enough to 
refer to the corresponding cigenmodc as the condensate 
wavefunction. The phase winding Q c f of the condensate 
wavefunction increases in steps of 2tt to a maximum of 6tt 
as the lattice is rotated faster and faster (see Fig. 112( a)). 
This corresponds to a maximum of L — 1 = 3 quantized 
vortices, with L x L being the size of the lattice. As 
in the single particle case, the maximum phase winding 
that can be observed within a lowest band model is lim- 
ited by the maximum phase difference of 7r/2 between 
neighbouring sites. 

Each vortex entry is associated with an avoided cross- 
ing between groundstate and first excited state. This is 
possible because both states have quasi-angular momen- 
tum, m = 0, allowing them to mix around the vortex 
entry point. This is demonstrated in Fig. 1 12( b) for the 
entry of the second vortex. 



B. Symmetry-incommensurate filling 

In behavior similar to that for four particles, the 
cigenmodc corresponding to the condensate wavefunction 
is macroscopically occupied, and Q c f for a symmetry- 



incommensurate number of particles increases in steps of 
27T up to a maximum of 2tt(L — 1). Yet, in contrast with 
a symmetry-commensurate filling, the discrete rotational 
symmetry of the system changes with each vortex entry. 
This is demonstrated in Fig. 113( a) and (b) where we plot 
both the phase winding and the quasi-angular momen- 
tum as a function of f£ for the case of five particles in a 
4x4 lattice. In this setting, the maximum phase winding 
is given by 6tt while the quasi-angular momentum takes 
the values m = 0, 1, 2, 3 in accordance with Eq. (|19fl . 

Since for symmetry-incommensurate filling the symme- 
try of the many-body wavefunction is different on cither 
side of the jump in the phase winding, transitions be- 
tween vortex states cannot occur via the mixing of energy 
cigenfunctions with the same symmetry. Hence, changes 
in vorticity are not associated with an avoided crossing 
between the ground state and the first excited state. In- 
stead, the transition occurs as the energy of an excited 
state with different quasi-angular momentum and phase 
winding drops low enough to become the new ground- 
state. The signature of vortex entry is thus a crossing of 
energy levels with different discrete rotational symmetry 
and phase winding. Fig. 113( c) depicts the level crossing 
associated with the entrance of a second vortex into a 
system of five particles in a 4 x 4 lattice. The level cross- 
ings are a non-trivial result for many particles since they 
correspond to a symmetry change in the ground state 
as a function of a parameter of the Hamiltonian and are 
indicative of quantum phase transitions |27| . 



VII. CONCLUSIONS 

We have studied zero temperature hard-core bosons 
in 2D rotating square lattices for filling factors of less 
than one atom per site using a modified Bosc-Hubbard 
Hamiltonian. An important feature of the system is the 
quasi-angular momentum, reflecting the discrete rota- 
tional symmetry of the lattice. Vortices enter the system 
as the angular velocity is ramped up. The number of vor- 
tices is obtained from the phase winding of the conden- 
sate wavefunction around the perimeter of the system. 
A lattice of size L x L can contain at most (L — 1) vor- 
tices in the lowest band model. We see quantum phase 
transitions as the quasi-angular momenta of the ground 
state changes. These are associated with vortices enter- 
ing a system which has filling incommensurate with the 
symmetry of the lattice. 

Even though we have studied small quantum systems, 
our work has implications in a broader context. The 
rotating lattice system is a promising experimental ap- 
proach allowing one to access more easily the regime of 
strong quantum-correlations, a major goal in the field 
in recent years. The novel aspect — the angular rotation 
frequency of an optical lattice — provides an additional 
parameter which in principle allows the experimentalist 
to explore a new axis of phase space. 

In order to extend our calculations to study systems 
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which arc of the size which will be more typical in exper- 
iments, it is necessary to go beyond the exact quantum 
ground state calculations. Although this is an important 
avenue for future work, the basic features of the quantum 
phase transitions and the emergence of the vortex lattice 
should persist in larger systems. Here, we have not ex- 
plored in depth the effects of high rotation which include 
the fragmentation of the condensate and the emergence 
of the physics of the fractional quantum Hall effect. In 
order to access this regime one should go beyond the low- 
est band description. 
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